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Abstract-This paper presents estimation of optimum speed 
parameters of wind turbine generator using new formulation 
of capacity factor. The formulation includes the effect 
coefficient of performance between cut-in and rated wind 
speed of the wind turbine generator. Second Mean Value 
Theorem for integrals and Weighted Average Harmonic Mean 
are applied to derive the capacity factor expression. The 
proposed model is tested using the actual wind data measured 
at Kappadgudda wind site. Analytically computed speed 
parameters are compared with the speed parameters of wind 
turbine generator installed at site to show the effect of 
coefficient of performance. It is observed from results that the 
accuracy in assessment of power in the wind and performance 
of wind turbine generator improves by 20 % with the new 
model. The proposed model helps in wind farm planning and 
site matching of wind turbine generator. 

Keywords- Coefficient of Performance; Capacity Factor; 
Wind Turbine; Power Curve 



speed ratio as shown in Fig. 1 m . In view of estimating 
higher energy production at higher capacity factor, this 
paper presents a derivation of an expression for CF 
considering the effect of aerodynamic properties of wind 
turbine blade indicated by C p . 

(^T^) Ideal C p without power regulation 
' 1 J Reduced C ^ by power regulation 




I. INTRODUCTION 

The average electrical power output depends upon wind 
velocities, turbine speed characteristics cut-in (V c ,) rated 
(V r ,) furling (V f ) and aerodynamic design of rotor blades. 
Capacity Factor (CF), tip speed ratio and Coefficient of 
Performance (Cp) are key factors that govern the optimum 
design of rotor blades and speed parameters of wind turbine 
generator (WTG). In order to estimate the performance of a 
WTG, appropriate modeling of CF considering aerodynamic 
analysis of blade and design parameter optimization for 
wind turbine is very vital. Aerodynamic design of wind 
turbine blades is dictated by Cp because Cp relates the wind 
speed, available mechanical power at the shaft and 
combines all the essential aerodynamic properties of a wind 
turbine. Further Cp is a turbine specific data and is normally 
provided by the manufacturer, during manufacturing of 
wind turbine, if accurate assessment of C p is not considered 
energy will be lost during the high wind velocities. 

Researchers have modeled CF to estimate average 
power production and performance evaluation of WTG 11-41 
assuming the value of C p is constant between V c to V r . The 
annual and monthly CF has been statistically computed 111 
with C p as a constant value. Modeling of CF in estimation of 
optimum speed parameters of wind turbine using 
normalized power curve 121 with assumption that Cp equal to 
Cp r . Further variation of Cp is taken in account while 
modeling CF r4 ~ 71 , but Cp is taken as Cp^ between V c and 
V r . However, Cp of wind turbine is neither constant nor 
maximum between V c and V r and it varies as function of tip 



O I ! ! , ^V(mA«J 

Fig. 1 Variation of Cp of WTG System 

H. PROPOSED MODEL 

The electrical power output of generator P e is given by 
the product of mechanical power output and the generator 
efficiency r] g 111 

where 

: air density (typically 1.225 kg/m 3 at sea level with standard 
conditions 

A : area swept by the rotor blades 
v : wind speed 

n m , rj g : mechanical and generator efficiency 
xoefficient of performance of the wind turbine is a function 
of tip speed ratio X and blade pitch angle /?. 

Aerodynamic efficiency C p of wind turbine is the ratio 
of turbine power to wind power 15 l The power coefficient C p 
is given in terms of the blade pitch angle /? and the tip-speed 
ratio X. The assumption of linear torque versus speed 
characteristics are not valid for wind turbines, which are not 
operated over wind range of wind speeds 15 '. Analytical 
analysis of aerodynamic system of such type of wind turbine 
is more complex. Consequently, numerical approximations 
have been developed in order to calculate the mechanical 
power characteristic of the wind turbine and a bi- 
dimensional characteristic function of Cp has been used and 
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validated in the laboratory . Further to simplify the 
approximation of C p , this paper considers blade pitch angle 
to be zero and is represented by non linear function as 



C , 



A, 



■ 0.244 



cor 



ri3o 



A 



6.56 \e 



(2) 



(3) 



v 



where co=rotational speed (rad/s), radius of rotor in m. 
CF is defined as the ratio of average electrical power 

[21 

output to the rated electrical power output 
1 



CF 



V, C, 



\C p V 3 f(y~)dv + ]f(y)d\ 



(4) 



, where C pr = Coefficient of performance at rated wind 
velocity, f(v) is the Weibull probability density function. 

Power extraction is higher, if variation of Cp is taken 
into account during modeling of CF as compared to the 
constant C p as shown in Fig. 2. 



Power curve using variation of C p 
Power curve using constant t p 




Wind Speed n n /s 
Fig. 2 P, versus wind Speed 

Substituting (2) in (4), the expression for CF can be written 
as: 



CF 



31.72 

V. 3 C.mr 



^e-" 3VUor v 4 f(v)dv 



(5) 



1.6064 



V 3 C r 



je- ,33V, ™'v 3 f(v)dv+jf(v)dv 

■ pr Vc V, 

Equation (5) contains three integrals and is solved using 
SM VT m . SMVT is applied to first integral of (5) to obtain 
extreme functions by substituting V r and V c in (6). 



CF 



corC (e 



13-3V c /o»- , „ — 13.3 V, 



31.72 



C pr (ar Vc 



-13.3V, /tor 



/(vj 



(6) 



-13.3V. I cor 



Weighted average for SMVT can be determined using 
Arithmetic Mean (AM), Geometric Mean (GM) and 
Harmonic Mean (HM), however in the paper HM is 
considered as 

• It is based on all observation 

• It is not affected very much by fluctuations of sampling. 

• It is rigidly defined 

• It is particularly useful in averaging special types of 
rates and ratios where time factor is varying. 



HM 



1e 



-\33(Vc+Vr)l cor 



(V) 



Substituting HM in first and second integral of (5) 
corresponding expressions are given by 



31.72HM 



b 



K C pr cor Vc 
1. 6064 HM 

v 3 c„ 



J v 4 f(v)d(v) 



v I 

J v 3 f{v)d(v) 



c = J/(v)rf(v) 



(8) 
(10) 

(11) 



Substituting (a), (b) and (c) in (5), modified expression 
of the CF is given. 



CF 



31.72HM 

V 3 C pr cor v 



jv 4 f(v)dv 



\.6Q6ArHM f 



V 3 C, 



jv 3 f(v)dv+ jf(v)dv 



(12) 

Further expression (12) is solved using integration by 
parts [2] detailed derivation is given in Appendix. New 
formulation for CF is given by. 



V, 



V. 



1 .6()64e-' 



— 13.3CV*-t-V<0 / 



c ,„- C«" 



, —1 3.3 V, 



, -(V r /c) k -(V Icy 

+ e r — e f 



(13) 



where V r /c: Normalized rated speed, T: Gamma function y: 
incomplete Gamma function. 



In Equation (13) CF is expressed in terms of 
aerodynamic properties of blade and normalized rated wind 
velocity. First two terms in the expression of CF contains 
weighted HM, which incorporates the variation of Cp into 
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the CF. Normalized average power and Turbine 
Performance Index (TP I) curve are defined as 



[21 



TPI 



V.nrV.r (PZ2-)AC 3 



P N *CF 
P * CF 

N ,max max 



CF 



(14) 



(15) 



CF and P N are computed using (13) and (14) for 
different values of V/c for known values of Weibull 
parameter. The normalized rated speed, V/c can be obtained 
corresponding to any point on the TPI curve. The value of 
rated speed V r for the site with known scale factor at TPI lmx 
becomes the optimum rated wind speed. Other two turbine 
parameters, namely cut-in and cut-out velocities, are 
estimated using the relations V c = pV r and Vf- qV r , where 
p<l and q>\. Typical values of p and q for Indian sites are 
found to be 0.275 and 1.85, respectively [2 l 

IE. RESULTS AND DISCUSSIONS 



Proposed CF model is validated using two case studies. 
SMVT is used to incorporate the variation of Cpinto the CF. 

Table I COMPARISON OF CF WITH AND WITHOUT Q, AT KAPPADGUDDA WIND SITE 



Monthly CFs are obtained from new formulation of CF 
using Cubic Mean Cube root (CMC) wind velocities. 
Normalized power, CF and TPI curves are drawn for the site 
whose Weibull shape k and scale c parameters known. 
Speed parameters for wind turbines are estimated and 
compared with installed turbine at site. 

Case Study-1: Kappadgudda wind power station, 
Kamataka state, India, is used for this case study. 
Analytically computed values of monthly and annual CF 
using proposed model, existing method [21 and actual 
capacity factors measured at Kappadgudda site are given in 
Table I and Fig. 3. The optimum turbine performance index 
curve, drawn for annual k= 2.3505 and c=9.8054 for wind 
site is shown in Fig. 5. The speed parameters derived from 
Fig. 4 are given in Table II. 

Case Study-2: Purpose of this case study is to validate 
the proposed model for the second set of wind data collected 
in the year 2007 from Kappadgudda wind site. Analytically 
obtained monthly and annual CFs with variations of C p , 
without Cp and actual CFs are given in Table I and Fig. 4. 
Speed parameters of wind turbine for the site are derived 
from Fig. 6, given in Table n. 



Data for the year 1998 


Data for the year 2007 


Month 


CMC 
m/s 


c P 


c 


k 


Capacity Factor 


CMC 
m/s 


c P 


c 


k 


Capacity Factor 


With 

c. 


Without 

Cp 


Actual 


With 

Cp 


Without 

Cp 


Actual 


Jan 


5.7955 


0.4006 


6.4796 


3.1085 


0.1067 


0.1073 


0.1377 


6.57 12 


0.4342 


7.309 


3.447 


0.158 


0.1680 


0.2307 


Feb 


6.3556 


0.4271 


7.1663 


2.4533 


0.2014 


0.1629 


0.1915 


6.7066 


0.4378 


7.511 


2.998 


0.186 


0.1900 


0.2414 


March 


6.0523 


0.4143 


6.8338 


2.2102 


0.1855 


0.1513 


0.1710 


5.8867 


0.4058 


6.644 


2.322 


0.143 


0.1502 


0.1973 


April 


7.8349 


0.4480 


8.8468 


2.1894 


0.3774 


0.2928 


0.1869 


7.4274 


0.4480 


8.386 


2.251 


0.296 


0.2801 


0.2254 


May 


10.9003 


0.3764 


122484 


2.7570 


0.6731 


0.5408 


0.4060 


9.9361 


0.4076 


11.121 


3.032 


0.582 


0.4932 


0.5284 


June 


10.4821 


0.3905 


1Ij6194 


3.6749 


0.6982 


0.5194 


0.5410 


10 5 936 


0.3868 


11940 


2.481 


0.620 


0.5361 


0.6543 


July 


11.7628 


0.3460 


13.1139 


3.2969 


0.7726 


0.6225 


0.5732 


1 1 .1 845 


0.3666 


12600 


2.533 


0.643 


0.5702 


0.6965 


Aug 


10.4108 


0.3928 


115997 


3.3358 


0.6736 


0.5107 


0.4592 


103815 


0.3938 


11570 


3.301 


0.637 


0.5434 


0.6900 


Sept 


8.2827 


0.4444 


9.3133 


2.7045 


0.4235 


0.3170 


0.3874 


9.3345 


0.4242 


10.400 


3.352 


0.523 


0.4461 


0.5510 


Oct 


5.6524 


0.3916 


6.3811 


2.2725 


0.1408 


0.1215 


0.2319 


6.6055 


0.4351 


7.455 


2.334 


0.215 


0.2012 


0.2123 


Nov 


7.5004 


0.4482 


8.3230 


3.6044 


0.2854 


0.2193 


0.1497 


6.8428 


0.4408 


7.363 


3.234 


0.189 


0.1945 


0.2035 


Dec 


6.6385 


0.4360 


7.44 14 


2.9325 


0.2043 


0.1673 


0.1858 


7.1612 


0.4458 


8.049 


2.741 


0.254 


0.2301 


0.2693 


Annual 


8.6892 


0.4097 


9.8054 


2.3505 


0.4643 


0.3601 


0.3018 


6.57 12 


0.4342 


7.309 


3.447 


0.473 


0.4152 


0.4485 
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Fig. 3 CFwith Covariation andconstant Cp for Vestas RRB at 
Kappadgudda 
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• Power estimated is 20% higher by wind turbine using 
proposed formulation as compared to existing method. 



NonmiMzpri Rated S])t?ert C\\ c) 



Fig 5: Turbine Performance Index for Kappadgudda site (1998 data) 
Table II Turbine SpeedParameters obtainedfrom Fig .4 and Fig.5 



Vestas-RRB turbine, D=30m 
(Wei bull parameters at site 
k =23505, c=9.8054, for 1998 data, 
ra ra ted=43 rad) 


EnerconE-30, D = 30m 
(Wei bull parameter at site 
k =2.469, c=10.170 for 2007 
data, io ra tai=43 rad) 


Para 


Variatioi 


Cons tan 


Actua 


Variatioi 


Cons tan 


Actual 


mete i 


of C p 


c„ 


of C p 


c„ 


Vc 


3.50 


3.2 


3.5 


3.41 


3.30 


2.50 


Vr 


12.75 


11.7 


13.5 


12.42 


11.39 


13.00 


Vf 


23.58 


21.7 


25.0t 


22.98 


21.01 


25.00 


CF 


0.47 


0.36 


0.30 


0.498 


0.415 


0.449 


Pn 


1.05 


0.80t 


0.941 


0.950 


0.780 


0.920 



Discussions: It is very interesting to note that, the 
analytical monthly and annual CF obtained from proposed 
model with Cp are higher than the existing model by 23 % 
and 17 % during 1998 and 2007 respectively. Further it is 
observed that from the first case study that, V,/c =1.30 at 
which PIC is maximum, CF and P N normalized power from 
proposed model is 0.42 and 1.05 respectively. These 
parameters obtained are higher than existing methods and 
actual wind turbine installed at the site. Similarly for the 
second case study, at PIC maximum, normalized rated speed 
is V/c =1.24 shown in Fig. 6. It is observed from Table II 
that, marginal variations of CF, P N and speed parameters 
between proposed method and actual wind turbine installed 
at site. However variation is more between proposed and 
existing method. Finally it is concluded that speed 
parameters estimated for wind turbine from proposed model 
will yield higher energy at higher capacity factor for both 
study years. 

IV. CONCLUSION 

Estimation of optimum speed parameter of wind turbine 
using new formulation of capacity factor is presented in the 
paper. New formulation of capacity factor is derived 
considering aerodynamic properties of wind turbine blade. 
Aerodynamic property of blade is governed by the C p of 
wind turbine. SMVT and weighted HM are applied to 
considering the variation of C p between V c and V r during 
modeling of CF. Estimated speed parameters from the 
proposed model are validated by comparing with existing 
method and speed parameters of Vestas RRB and Enercon 
E-30 wind turbine installed at the site. Conclusions derived 
fromthe two illustrative case studies are as follows. 
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Fig 6: TurbinePerformance Index for Kappadgudda site (2007 data) 

• Speed parameters estimated for wind turbine taking into 
account the variation of C p will yield higher energy at 
higher capacity factor. 

• New formulation of capacity factor is a useful tool in site 
matching of wind turbine generators for wind potential 
site. 
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APPENDIX 

DERIVATION OF NEW FORMULATION OF 
CAPACITY FACTOR (CF) FOR WIND TURBINE 



A.l 



[5] 



Capacity factor of WTG is given by l2] 

, v. v f 
CF=— j Cp v 3 f(v)dv + j f(y)dv 

Coefficient of Performance (C p ) for WT is given by 

( 130 ^ -133/ 

C p = 0.2441 ^ - 6.56 \e /x 

Substituting (A.2) in (A.l), CF is obtained 

CF = ] 1J2 ] e- 133v ' " V/( v)dv - 1 4° 6 ^ I e- ,33V 1 "V/(v)rfv 



A.2 



■ \f{v)dv 



A3 

Second Mean Value Theorem (SM VT) for integral theorem 
states that 

b 

Y = j f (x)<I>(x)dx 

a 

A.4 

b 

Y = f(c)^<S>(x)dx 
A.5 

fic) is weighted Harmonic Mean (HM) r9 ' and is obtained 
by substituting upper and lower limits of (A.5) 

2f(a)f(b) 



f(a) + f(b) 



/(c) = HM 

A.6 

SVMTis applied to (A. 3) and weighted HM is given by 

2 e -13.3(Vc+Vr)/ojr 



HM = 



A.l 



e ~\3.3Vc I cor + e ~\3.3Vr/ cor 



Substituting (A.7) in (A.3), CF is given by. 



CF = 



v,/^ -13.3V, /-r, J v V'(v)dv 



„V c -13.3V. 

K C pr cor(e 

3212e -,3.3,V, + V^, r Vr V 

~ V 3 -,3.3V,. Im .3.3V,/ J V f^ V+ J 

v r l ">rV e ~ r e > Vc V r 

(A. 8) contains three integral 



63.44e 



-13.3(V r +V c )/ar 



' V r 3 C pr cor 



.we ' - f 4 
v f(v)dv, 

( -11,3V c I car _-13.3F, / (Or n J - 1 V ' 



A.9 



3 2i2e~ 13 ' 3(V '' +v, ' )/ ' or 



' != v,' C ; ( 7-'»-" + .-*<-> t' mdv - 

A.10 



■ 

7 3 = \f{v)dv 

v r 

A.ll 

Substituting Weibull probability density function in (A.9- 
A.ll) 

1 v/c^e- 1 ^ / -+^' / ' ,r )| v cUJ e V 

A. 12 

Using integration by substitution to (A. 12) l2] 

iyj ' then dx= ^/ki v /J~ 1 dv and 



X = < 



1/k 



V = CX 

A.13 



where: 



a = - 



63A4e~ 133(K+K)/mr 



C vr cor(e 



-1 3.3V C / cor + e -\33V r lcor y 



A. 14 



V,. 



Applying integration by parts (A. 14) 



63.44V r e-' 33iV ' +v -"" r 
corC(e' i " V '' r "+e- i33V '' r ") 



4r — | 

p * e -<y,io* _ e - ( v r /c,« + A^, r 
ft 



VX 4 
cj'k 



A. 15 



Similarly final expression is obtained for/ 2 and I 3 



l.(,Q64e- {3MVnVc, " ar 
C^e™' 1 " +e- l33v '"" r ) 



/ 3 =e -^/o' 

Combining (A. 15), (A. 16) and (A. 17), we get 
63.44V e - nxv ' +v < )lm 



A.16 



A. 17 



CF 



G)rC(e~ U - 3V ' lar +e~ n " V ' lm ) 
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1 .6064e 



-13.3(Vr+Vc)/(Dr 



-13.3V,. I cor 



-13.3V, I cor ^ 

e ) 



3r 



v,4\ r \[c J ■ u 



.g-C^/c) 4 _ c -(V / /c)« 



A.18 
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